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Abstract 

WMAP5 and related data have greatly restricted the range of acceptable 
cosmologies, by providing precise likelihood ellypses on the the w -w a plane. 
We discuss first how such ellypses can be numerically rebuilt, and present 
then a map of constant-^ models whose spectra, at various redshift, are 
expected to coincide with acceptable models within ~ 1 %. 
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1. Introduction 

One of the main puzzles of cosmology is why a model as ACDM, implying 
so many conceptual problems, apparently fits all linear data in such unrivalled 
fashion [10, H|. 

It is then important that the fine tuning paradox of ACDM is eased in 
cosmologies where Dark Energy (DE) is a self-interacting scalar field <fi (dDE 
cosmologies), with no likelihood downgrade 0, 0|. Unfortunately, however, 
only cosmological observables can provide information on the form of the self- 
interaction potential, even though several researchers incline to priviledge 
potentials allowing tracking solutions. 

When aiming to obtain information on the physical potential, the basic 
observable is however the evolution of the DE scale parameter, to (a). Here 
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a = 1/(1 + z) is the scale factor of a spatially flat metric 

ds 2 = c 2 dt 2 -a 2 (t)d£ 2 , (1) 

while z is the redshift. 

The analysis of available data, made by the WMAP team jjjj, was able 
to constrain the coefficients wq and w a in the expression 

w(a) = w + (1 - a) w a , (2) 

putting again into evidence that a model with w = — 1 is close to top likeli- 
hood, but also stressing a preference for the phantom area (w < —1), which 
is hardly consistent with current tracking potentials. 

The main tool, to go beyond these constraints on the w(a) law, will be 
tomographic shear analyses [6|, able to reconstruct the spectrum of density 
fluctuations at various z's, with a precision approaching 1 % |7|. 

This work is therefore focused on the relevance of spectral predictions, and 
aims at providing a tool to ease the determination of model spectra. Within 
this context, it is essential to outline the spectral equivalence criterion (SEC). 

It has been noted since a few years [g, ^ that the density fluctuation 
spectra P(k) = (\5p/p\ 2 ), up to k = 3/iMpc -1 , essentially depend on the 
distance from the LSB (Last Scattering Band). This was first verified at 
z = [lo| and then tested at higher z 11\. iV-body simulations show that 
the SEC is however true when model parameters are suitably tuned. 

Let us be more specific on this point: When a given w(a) is considered, 
it is easy to determine the comoving distance of the LSB. Keeping then the 
same values of f2& c and h, we can seek a constant-u> cosmology with the same 



distance from the LSB; [10J tested that the spectra of these two cosmologies 
coincide (within 1 %). They also saw that, when 2^0 values are explored, 
spectral discrepancies are mostly greater, in the range of a few percents. 

The SEC however works also at z ^ 0, provided that the distance between 
such z and the LSB is evaluated and a constant-^ auxiliary model is defined, 
with an equal distance between z and the LSB. The assigned cosmology and 
the auxiliary model must also have equal f2f, jC , h at z = 0, while <J%(z = 0) 
must be tuned in order that, at z, the r.m.s. density fluctuations of the two 



models, on the 8 /i _1 Mpc scale, coincide ll| . 

Let us outline, in particular, that the SEC does not require that the 
auxiliary model shares the values of c and h at the assigned z ^ 0. In fact, 
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Figure 1: 1- and 2-a curves, yielding the marginalized model likelihood on the WQ-w a 
plane, as obtained from the Lambda NASA site. The reproduction device is detailed in 
the text. Also 0.5- and 1.5-er curves are provided. Crosses indicates models for which the 
SEC was explicitly tested at various z's. 



by multiplying the critical density definition by fl C) b, one has 

u) c>b oc Q c , b H 2 = (87rG/3)tt Ctb p cr = (87rC7/3)p Cj6 oc a" 3 



(3) 



so that the assigned cosmology and the auxiliary model share the values of 
u C)b at any z, and this is enough. 

This comes as no surprise, however: most linear feature, e.g. BAO's, 
essentially depend just on uj b>c . 

This recipe, however, prescribes a different constant--?/; at any z and is 
somehow curious that, starting from the assigned w(a), one builds a w e ff(a) 
law, completely different from it. In a similar way, one can build a cr 8 .( 2=0 ) 
dependence from z, in order that auxiliary models be fairly normalized at 



any z. Examples of w e ff(z) and cr 8 ^ z=0 )(z) are provided by [U 



The scope of this work, instead, amounts to applying the SEC to models 
with w(a) of the form ([2]) and consistent with data, exploring decreasing 
likelihood ellypsoids. 

The plan of the paper is as follows: In the next Section we shall discuss 
how one can rebuild the likelihood ellypsoids on the WQ-w a plane. This 
Section defines the parameter t used in the plots which are one of the results 
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of this work. In Section 3 we shall then discuss such plots. Finally, Section 
4 will be devoted to drawing our conclusions. 

2. How to reproduce likelihood curves with Bezier paths 

Solid lines in Figure [1] reproduce the likelihood ellypsoids on the WQ-w a 
plane of Figure 14 in It is significant to explain how such reproduction 
is obtained, as the technique used is also essential to build the forthcoming 
Figures, where a parameter t appears, which are one of the results of this 
work. In a sense, in this way we pass from an eulerian w -w a description to 
a lagrangian t description, also stressing the simmetries of the likelihood. 

Dashed lines in Figure [1] approximately yield 0.5- and 1.5-er's; they are 
defined as the locus of points of equal t exactly at the center of the intervals 
between top likelihoods and 1 a or 1- and 2-a's. 

The technique we shall describe to draw curves on a plane, named after 
Bezier, is largely used in vector graphics to model smooth curves which can 
be scaled indefinitely, without any bound, by the limits of rasterized images. 

For instance, imaging systems like PostScript, use cubic Bezier curves: 

B(u) = (1 -u) 3 P + 3(l -u) 2 uP 1 + 3(1 - u)u 2 P 2 + u 3 P 3 we [0,1]. (4) 

The vector B runs on the WQ-w a plane, describing a curve when u varies from 
to 1. The curve is fixed by the positions of the points (k=0,...,3). In 
eq. flJJ), Pfc indicates a vector ending on the very P& point. The curve starts 
at P Q and is initially directed toward Pi] however, it bends soon, owing to the 
setting of all other points; for u — 1, it meets P3, its final direction being set 
by Pi. Clearly, it hits neither P\ nor P2. these points only provide directional 
information, while the distance |Po-Pi| tells us how persistently the curve 
moves towards Pi; an analogous effect, at the end of the run, is fixed by the 
distance |P 3 -P 2 |. 

Quadratic and cubic Bezier curves are mostly used; when more complex 
shapes are needed, rather than making recourse to higher degree curves, 
numerically expensive to evaluate, lower degree Bezier curves are patched 
together. In our specific case, the ellypsoid is shared in four paths, each 
described by an expression (j4]). Let us then label the paths with i; we can 
follow the whole ellypsoid with a single parameter t e [0, 1] by labelling the 
paths with an index i (i = 1, .., 4) and setting u — At — i + 1; accordingly, it 
shall be 

i — 1 < At < i in the i — th path. (5) 
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Using the PostScript file of Figure 14 in [3j, we can easily obtain the x, y 
coordinates of the Pk points of each path. In table 1 we report such points 
for l-o" and 2-cr ellypsoids, converted into WQ-w a units. 
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-1.3593 


1.2597 


-1.3765 


1.1038 


-1.2247 


0.4538 


-1.1387 


0.1068 
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-1.1387 


0.1068 


-1.0528 
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-0.9333 


-0.6217 


-0.9145 


-0.4527 
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-0.9145 


-0.4527 


-0.8942 


-0.2704 


-0.9542 


0.1407 


-1.0326 


0.4713 
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-1.0326 


0.4713 


-1.1318 


0.8893 


-1.3384 


1.4488 


-1.3593 


1.2597 
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-1.4967 


1.6957 


-1.5229 


1.5197 


-1.3801 


0.8544 


-1.1351 


-0.2619 
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-1.1351 


-0.2619 


-0.9558 


-1.0784 


-0.7862 


-1.6310 


-0.7469 


-1.3819 
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-0.7469 


-1.3819 


-0.7039 


-1.1091 


-0.8015 


-0.2921 


-0.9625 


0.4120 


4 


-0.9625 


0.4120 


-1.1235 


1.1160 


-1.4676 


1.8905 


-1.4967 


1.6957 



Table 1: Points denning the 4 cubic Bezier expressions yielding the 1— a and 2-cr curves 
(upper and lower table, respectively). 

These coordinates were used to produce Figure [TJ This description is 
however essential also to introduce the parameter t, which shall be used in 
the forthcoming plots. 

3. Results 

In fact, by letting the parameter t to run, we move along each curve of 
Figured], and also define the so-called 0.5-cr and 1.5-er curves. 

In Figure [2] we report the variations of w (dashed curve) and w a (dotted 
curve), when t runs from to 1 (we plot w a /3 — 1 instead of w a , to allow 
wider ordinate spacing). Each t value, then, yields a WQ-w a couple, defining 
a model at 0.5-cr from top likelihood. In the same Figure, we also report 
how the constant state parameter w e ff depends on t, at various redshift, 
for a model yielding the same distance of 0.5-cr models from the LSB. The 
redshift values considered are z = 0, 0.5, 1, and 2 . When z increases, the 
ordinate interval spanned by w e ff values becomes wider; this is true at any 
number of cr's and allows to individuate the solid curve referring to each z. 
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Figure 2: This plots gives the state parameters for constant-w models (solid lines), yielding 
the same distance from the LSB of Wo~w a models at 0.5 ct's from the top likelihood 
cosmology; the 4 solid line refer to z = 0, 0.5, 1, and 2 . 
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Figure 3: As Figure for l-er models. 
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Figure 4: As Figure^ for 1.5-cr models. 
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Figure 5: As Figure [2J for 2-er models. 
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Figures [3H5] provide analogous results for models at 1, 1.5 ad 2 a from 
the top-likelihood cosmology. 

These figure are obtained by taking uj c = 0.228, cj& = 0.046, h = 0.71. 

Let us notice that solid curves converge and meet dashed ones for t values 
where w a vanishes. These points also set a transition between the t-intervals 
where w e // increases or decreases with z. 

The top likelihood, as is known, corresponds to a phantom-DE model. 
At 0.5 cr's from it most models have wq < —1 but, also a large fraction of 
the few models with u> > —1 correspond to w e ff < —1. On the contrary, 
there are quite a few models, characterized by wo < —1 whose spectra are 
equivalent to constant-it; models with w e jf > — 1. 

The t interval where w$ < — 1 becomes wider when a greater number of 
cr's is considered. On the contrary, the width of the intervals characterized 
by wo < —I and w e ff > — 1 does not change much with the number of a's. 

4. Conclusions 

When tomographic cosmic shear data will become available, an inspection 
on DE nature will surely start from comparing them with the predictions of 
constant-if cosmologies. 

As soon as data will become more refined, it will be possible to bin them, 
discriminating among the w values best fitting data at various redshift. 

It is quite possible that these inspections yield w values compatible with 
a constant, all through the redshift range explored. As it is possible that 
such value is compatible with —1, so vanifying the efforts to improve our 
knowledge on DE nature. 

Let us suppose that, instead, data analysis is consistent with the same 
Ub tC values in all bins, but require different w 's in different bins. We wish to 
outline here a major danger that data analysis could meet in this welcome 
case. As a matter of fact, one could be tempted to conclude that the w(z) 
dependence found is the physical scale dependence of the DE state parameter. 
Unfortunately, this could be badly untrue. 

In fact, cosmic shear spectra can be easily translated into fluctuation 
spectra, so that observational values of w(z) would correspond to w e ff(z), 
not to the physical w(z). 

How different the two behaviors can be is already represented by Figures 
2-5 in this paper. But Figures [6] and [7] further illustrate this point. They 
show what is w e ff when w has a given value, for models laying on the 1—a or 
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Figure 6: Bias in observational w values (named here w e f /), in respect to the physical value 
of the state parameter, at each redshift. The points in the curves of the upper frame refer 
to different models at 1-a from the top-likelihood cosmology, discriminated by different t 
values. In the lower frame the ratio w e ff/w is explicitly shown as a function of t . 
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Figure 7: As Figure [6] for models at 2-a from the top-likelihood cosmology. 
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2— a curves. Different colors in each plot correspond to different z's. In the 
lower frame, the w e ff/w ratio is also shown, as a function of t. Even in the 
most furtunate cases, provided that the true cosmology is not a constant-^ 
one, discrepancies are hardly below 10 % and are greatest at z approaching 
zero. 

Another point concerning model fitting is that a systematic mapping 
of (reasonable) dDE cosmologies, in order to fit future data, is apparently 
unnecessary. At each z, in fact, there will be a constant-w model able to fit 
any dDE cosmology. dDE will then be revealed by the w dependence on z, 
as above outlined. At present, well approximated analytical expressions of 
P(k), at different z, are provided by the so-called Halofit formulae, holding 



for ACDM cosmologies [12J. Some attempt to generalize Halofit to constant 



w ^ — 1 were also performed [13j , but they do not cover the desired parameter 
ranges. Our conclusion is that it will be important to extend Halofit to 
constant-w; cosmologies, for the whole range of (reasonable) cosmological 
parameters; this will enable us to fit future cosmic shear data without any 
substantial restriction on the w(a) behavior. 
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